Understanding B —> J / in the Standard Model 



Ying Li* 

Department of Physics, Yantai University, Yantai 264005, China 

Hai-Yang Cheng 

Institute of Physics, Academia Sinica, Taipei, Taiwan 115, R.O.C. 

May 20, 2009 



Abstract 

The rare decay B d — ► J/y<p can proceed via four distinct mechanisms: (i) production of 

the via tri-gluon fusion, (ii) photoproduction of the J/y or 0, (iii) final-state rescattering of 

(*) (*) 

D s D s produced in the B c \ decay to J/yty, and (iv) production of the via ft) — (f> mixing. 
In this work, we examined the contributions of photoproduction and final-state rescattering to 
B° d — > J/y<j) and found that the corresponding branching ratios were of the orders 10~ n and 
10~ 9 , respectively. Hence, this decay is dominated by the ft) — mixing effect. 

1. The observation of B decays to charmonium provides important evidence for the Cabbio- 
Kabayashi-Maskawa model, as well as an important advance in our understanding of the Standard 
Model and QCD dynamics. Recently, Belle reported an upper limit 9.4 x for the branching 
ratio of B° — > J/\j/(j) at the 90% confidence level JT). This process is expected to be suppressed by 
the Okubo-Zweig-Iizuka (OZI) rule [2] disfavoring disconnected quark diagrams. 

The main processes for 5° ^J/y<$> can be sorted into four different classes: (i) the neutral 
vector meson is produced through tri-gluon fusion (Fig. [Q), which is formally the reason why 
this channel is OZI-suppressed, (ii) the J/\y or arises from a photon emission, followed by 
fragmentation (Fig. ©, (iii) the decay particles J/xy and are produced through long-distance 
final-state interactions (FSI) (see Fig. [3]), and (iv) the § comes from the decay of B — > J / i//<y 
followed by ft) — ^ mixing; that is, is not a pure ss state and contains a tiny qq component. 

In ||3) , Gronau and Rosner pointed out that the major contribution to the decay B d — >• J / xyty 
arises from ft) — § mixing. Neglecting isospin violation and the admixture with the p° meson, one 
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mechanism. 



can parameterize mixing in terms of an angle 8 such that the physical 0) and are related to 
the ideally mixed states co 1 = (uil + dd) / y/2 and § l = ss by 

co \ I cos 8 sin 5 \ / (0 1 
(j) J I — sin 8 cos 8 I \ (j) 1 

and the mixing angle is approximately 8 = —(3. 34 ±0.17)° H|. Within this mechanism, the 
authors estimated the rates of this decay mode and other similar processes in 5° and B® decays, 
and found that the Belle's upper limit is about a factor of five above their estimation. Also, they 
argued that the final-state rescattering contributions to this decay mode are very small and can be 
neglected. 

Let us make crude estimates of the various contributions to B — > J/yf(j> by the aforementioned 
four mechanisms. Due to the complicated QCD dynamics, it is difficult to calculate the tri-gluon 
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fusion reliably. Roughly, the tri-gluon fusion contribution gives 



&(B° -> J I yQ) tri _ gluon = ^(B° -> J/\jfa>)a* « 2.7 x 1(T 5 x (0.3) 3 ~ 2.4 x 1(T 8 , (2) 

where use of ^(5° -> « -> //yp) = (2.7 ±0.4) x 10 -5 [5] has been made. The 

contribution of photoproduction is calculable to the leading power of the 1/m^ expansion and is of 



^(B° - //V^photoproduction = ^(B° - J/vy)o£m ~ 1(T 7 x (1/137) 2 ~ KT 11 . (3) 



In the final-state rescattering picture, the B — * //V0 decay proceeds via a 5 meson decay into 
D s D s through W-exchange followed by a rescattering of D s As to J/ through D s 
exchange. It is anticipated that 



£g(B° -> //v/0)fsi = =^(5° -> Dr )+ Dr )_ ) (lCT 3 - 1(T 4 ) ~3x 1(T 8 - 3 x 1(T 9 , (4) 



where the analysis of final-state interactions in B — > (j)K*,pK* suggests that the rate of the B- 
meson decay into the final state under consideration (for example, 5° — > 7/ty/0) is suppressed 
relative to that of the intermediate state (B° — > in this example) by three to four orders 

of magnitude 0. Finally, the production of// i//0 through co — mixing is expected to be 

^(5° -> 7/^0)(b-0 mixing = ^(5° -> sin 2 5 « 2.7 x 1(T 5 x (0.08) 2 ~ 1.7 x 1(T 7 . (5) 

Therefore, the rare decay B — > J /yty is indeed dominated by the <y — mixing effect. 

In this letter, we will study the effects of photoproduction and final-state rescattering in more 
detail even though they are not the main contributions to B — > J /yty. We wish to have quantitative 
results to confirm the above crude estimates. 

2. Firtst, Let us evaluate the photoproduction, which plays an important role in decay modes 
such as B — > p^*, p0 0. In this mechanism, B®, — > J /yty can be regarded as the cascade process 
B Q d —> J /yy ^ J /yty or 5° — > 0y — > J/yr^). The radiative decay 5 — > Vy has been well studied 
in the frameworks of the QCD factorization approach [8], the perturbative QCD approach (pQCD) 
[|9l and soft-collinear effective theory ifTOl . Due to the suppression of Wilson coefficients, we will 
neglect the contribution from B~ d — > 0y — > J / 

According to the Feymann diagrams depicted in Fig. [2l the amplitude of B~ d — > J/tyty can be 
written as 



order 




(6) 
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where we have used (0\sY v s\^)) = — m^/^£ v * and and are the decay constant and mass of 
the meson, respectively. Therefore, we obtain the result 



R 



3m„ 



0.0003, (7) 



with fa = 0.237 GeV. In the literature, it has been estimated that &{B° d -> //V7) = 7.7 x 10~ 9 
ifTTIl in QCD factorization and B8{EP d -> 7/y/y) = 4.5 x 10~ 7 [12] in perturbative QCD. Therefore, 
the predictions of QCDF and pQCD differ by one to two orders of magnitude. The possible reason 
for this huge discrepancy was explained in Ref. lfT2l . Roughly speaking, this is due mainly to the 
use of different 7/i/a wave functions in Ref. lTTTII and Ref. lfT2ll . If the charm quark is heavy, the 
wave function of J/xy will be symmetric under x <-» 1 — x and sharply peaked around x = 0.5. 
However, the cross section of e + e~ — > T] c +J/\ff calculated within the NRQCD approach is much 
smaller than the experimental data. Bondar and Chernyak lfT3l have pointed out that the origin of 
the discrepancy is due to the fact that the charm quark is not heavy enough and, as a result, the 
charmonium wave functions are not sufficiently narrow for a reasonable application of NRQCD 
to the description of charmonium production. Using more realistic models, these authors have 
proposed a new wave function for 7/ Xjf, which can be used to explain the data well. This new wave 
function is employed in Ref. [12], while the delta function is used in Ref. [1 1]. 

Even taking the pQCD result for B d — > J/xyy, the photoproduction mechanism leads to a very 
small branching ratio for B~ d — > J/xj/Q of order 10~ H , which is not accessible even at the future 
Super-B factories. Since the is produced from a virtual photon which is transversely polarized 
mostly, the longitudinal polarization of the decay B — > 7/v0 via photoproduction will be very 
small. 

3. As mentioned above, 5° — > J '/xj/Q receives long-distance contributions from a B meson 
decay into D ( * )+ D { * ] ~ followed by a rescattering of D { * )+ D { * y to J/\ff0. The D^D^ states 
from B decays can rescatter to J/yrfy through the f -channel D; s exchange in the triangle dia- 
grams depicted in Fig. [3] Before proceeding, we would like to remark briefly on the motivation for 
considering the rescattering mechanism with Dl exchange. At the hadron level, final-state inter- 
actions manifest as the rescattering processes with s-channel resonances and one particle exchange 
in the f-channel. Due to the lack of the existence of resonances at energies close to the B meson 
mass, we will therefore model FSIs as rescattering processes of some intermediate two-body state 
with one particle exchange in the /-channel, we will compute the absorptive part via the optical 
theorem [6]. We consider charm intermediate states based on the idea that if the intermediate states 
are CKM more favored than the final state, then the absorptive part of the final-state rescattering 
amplitude can easily give rise to large strong phases and make significant contributions to the rates. 
It has been shown in Ref. [[6] that the direct CP- violating partial rate asymmetries in charmless B 
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decays to kk/kK and pn are significantly affected by final-state rescattering and their signs are 
generally different from those predicted by the short-distance approach. Especially, the calculated 
CP asymmetry Acp(K + n~) = — 0. 141qq3 for B° — > K + n via rescattering [[61 agrees with experi- 
ments in both magnitude and sign, whereas the QCD factorization prediction Acp(K + 7t~) ~ 0.045 
lfl4l is wrong in sign. This example illustrates that the rescattering approach gives a reasonable 
description of FSIs. 

To evaluate Fig. [3j we note that the effective Lagrangian for QDyDy vertices can be found 
in J6l, and the effective Lagrangian for J / y/D^ vertices is given by 

X w d s d s = ig V D s D s ¥n{d^D s D;-D s d^D s ^, (8) 

^ v d*d s = -2f V D*D s e^d^ V (daD s *pDj+D s d a D*^, (9) 

+xif v DTd^ - xir v d^ v D s *^}. (10) 

(*) (*) 

The coupling constants for the (j)D s D s vertices can be related to the parameters gy, j5 and A 
appearing in the effective chiral Lagrangian describing the interactions of heavy mesons with low 
momentum vector mesons [|20l in the following manner 



8<j>D s D s = ^j=- = 3.75 , f<j,D s D* s = ^j=- = 2.30 GeV l , 

A 

f<j>D*D* = ~^ mD " s =4-85, g w d sDs = 4fy D * Ds = 4fy D * D * /m D * = 10, (11) 

where we have assumed /3 = 0.9 and A = 0.56 GeV -1 ETTl and the relation gy = m p J f n EOl . The 

couplings for J/\\td[ *'D$ are taken from Ref. ll22l based on an effective field theory of quarks 

(*) (*) 

and mesons. Note that the same (j)D s 'D s vertex also appears in the rescattering contribution to 
B — > (f)K*. A study in [6J shows that the rescattering mechanism via D s exchange can enhance 
the rate and yield a large transverse polarization in B — > (j)K*. 

o (*)— (*)+ 

In total, there are eight different FSI diagrams in Fig. [3J The B — > D s D s — > J/y/Q 
amplitudes via Dg exchange are similar to the B — > D^D^ — > K*§ amplitudes via Ds exchange 
that have been studied in Ref.[6]. Therefore, the amplitudes of the former can be obtained from 
the latter through the replacements K* — > J/\]f and dW — > D s *\ For example, the absorptive part 
contributions of B° — > DjDf —> J/ yf(j> amplitudes via D s exchange is given by 

A b s(D 7 Di;D,) = 1/ ^ ^ (2. W,. - » - P2 )A(? D 7 Dl) 

x (2i)g DAt F(f "' k)F i P2 ' k) {-2i)g DM/v (%* ■ Pi) W • Pi), (12) 

' m D s 

where k = p\ — p^ = P4 — P2 is the momentum of the exchanged particle. Since the particle ex- 
changed in the t channel is off shell and since final state particles are hard, form factors or cutoffs 
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Figure 3: Long-distance contribution to B — > J/y/cj) 



must be introduced to the strong vertices to render the calculation meaningful in perturbation the- 
ory. The form factor F(p,k) for the off- shell effect of the exchanged particle can be parametrized 
as 



normalized to unity at t = ml xc , where m exc is the mass of the exchanged particle. The cutoff A 
in the form factor F(t) should be not far from the physical mass of the exchanged particle. To be 
specific, we write [0 



where the parameter r\ is expected to be of order unity and it depends not only on the exchanged 
particle but also on the external particles involved in the strong-interaction vertex. As we do not 
have first-principles calculations for form factors, we shall use the measured decay rates to fix the 
unknown cutoff parameters. Although the strong couplings are large in magnitude, the rescattering 
amplitude is suppressed by a factor of F 2 (t) ~ (m Aq CD /Y ) n . Consequently, the off-shell effect 
will render the perturbative calculation meaningful. It is also evident from Eq. (fT2l) that the final- 
state rescattering contributions vanish in the heavy quark limit, as it should be. 

As discussed in Ref.[6], the FSI contribution from the B — > DjDf decay will affect both At 
and An amplitudes of the 5 — > // decay, whereas both B — > D*D S and B — > D S D* will affect only 
the A_l term of the B — > J/yrQ decay amplitude. Finally, the FSI effect from the decay B — > £>*£>* 
contributes to all three polarization components A/im. 

In order to perform a numerical study of the long-distance contributions, we need to specify the 

o (*)+ 

short-distance A (B —> D s As ) amplitudes. This decay proceeds only through VK-exchange, 
and it can be calculated in pQCD effectively without introducing any new parameters lfT51 [T6ll . 
Numerically, we have (in units of V c bV* d GeV) 




(13) 



A = /n e xc + T?AQCD, 



(14) 



stf{B»^DtDJ) 



7.93 x l(T 6 + z'0.94x 10 



-6 



0.98 x 10~ 6 + zT.12x 10~ 7 , 



(15) 
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and 



a 



1.9 x KT 6 -?'1.4x 10" 7 , 



b = -6.5 x 10~ 9 + *4.7x 10~ 8 , 

c = 6.7 x l(T 9 -?4.9x 10~ 8 , (16) 

for the 5° — > amplitude given by 

^(B -> D;+( Pl!ei )Z)r (p 2 ,e 2 )) = a(e? • ef) + fc(e? ■ p 2 )(e| -pi) + /ce a ^ v e 1 *"e 2 *^M- ( 17 ) 

It follows that the branching ratios of B — * ; D y s ' read 

^(B°^D+D-) = (3.3 ±1.1) x 1(T 5 , 
@(B°->D*+DJ) = (2.6 ± 1.0) x 10~ 5 , 

gg(B° D*+D*-) = (1.2±0.4) x 10~ 5 . (18) 

In the above calculation, we have included the errors coming from the the hadronic wave functions 

(*) 

that are dominated by the D s meson distribution amplitude rather than the B meson, as the latter 

is more or less fixed by the well measured channels such as B — >• K%, %%. Since we employ the 

u\ 

updated D s distribution amplitude 11171 
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(x, b) = —f nM x{ \-x) 1 + A nW ( 1 - 2x) 



V6 Jd 



Di 



with a (») = 0.5 GeV and (O = (0.6 — 0.8) GeV, our predictions are slightly smaller than the ones 
in fTTSll H but consistent with the current experimental limits |[T8l[T9l 

&(B° -> D+DJ) < 1.0 x 10~ 4 (BaBar), < 3.6 x 10~ 5 (Belle), 
^(B° D* + Dj) < 1.3 x 10" 4 (BaBar), 

^(B° -> D* + D*-) < 2.4 x 10~ 4 (BaBar). (20) 

For the parameter 77 in Eq. (fl4|) . we shall use the one T] =0.80 extracted from B — > (j)K* [|6]|. 
With the B ^ Dl ' D^' amplitudes given before and the parameters (Ull) . the decay rate and the 
longitudinal polarization fraction fa of B — »• 7/ due to final- state rescattering turn out to be 

38(eP -*J/y$)?si = (3.7^) x 10 -9 , / L = 0.41±0.02. (21) 



'Our estimate of £$(B° — > DfDj ) is smaller by more than a factor of two than a value of (7.8+j'g) x 10~ 5 obtained 
in lfl5l using the same PQCD approach. This is mainly due to the additional exponential term exp(— co 2 b 2 /2) in the 
revised Dg distribution amplitude, Eq. ( fT~9b . Based on the diagrammatic approach, an estimate of £$(B° — > DfD~) = 
(4.01} j) x 10- 6 was obtained in Ref. 1231 . which is smaller than the PQCD result by one order of magnitude. This 
should be checked by experiment. 
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Here we only show the major errors stemming from the uncertainties in the parameter r\ and the 
cutoff scale A (see Eq. (fl4l) ) where we have assigned a 15% error to Aqcd an d an error of 0.01 
to 77. As in Ref.[6J, we have assumed monopole behavior [n = 1 in Eq. (fT3l) l for the form factor 
F(t,mo s ) and a dipole form (n = 2) for F(t,mD*)- It should be stressed that the estimate of the 
FSI contributions is model-dependent as it depends on how we model the final-state rescattering. 
In view of this point and the theoretical discrepancy between PQCD and the topological diagram 
approach for the rate of B — > D s D s , it is conceivable that the actual theoretical uncertainties 
are considerably larger than those given in Eq. (T2T)) . At any rate, it is evident that the final- 
state rescattering contribution to — > J /yty is smaller than the effects of 0) — mixing by two 
orders of magnitude. We thus confirm the argument by Gronau and Rosner [3] that a significant 
enhancement of this mode by rescattering is unlikely. 

4. In this work we have examined the contributions from photoproduction and final- state 
rescattering to B d — > J/\j/(j> and found that the corresponding branching ratios are of order 10~ n 
and 10~ 9 , respectively. Hence, this decay is dominated by the 00 — mixing effect as advocated 
by Gronau and Rosner. 
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